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Abstract 



Previous calculations of the polarizabilities of the nucleon within 
the framework of heavy baryon chiral perturbation theory have in- 
cluded the contribution of the A(1232) only in its effects on various 
contact terms or have been performed in chiral SU(3) where system- 
atic errors are difficult to control. Herein we perform a corresponding 
calculation in chiral SU(2) wherein A(1232) is treated as an explicit 
degree of freedom and the expansion is taken to third order in soft mo- 
menta, the pion mass and the quantity M^ — Mj^, collectively denoted 
e. We present the results of a systematic 0{e'^) calculation of forward 
Compton scattering off the nucleon, extract the electric polarizability 
(Xe, the magnetic polarizability Pm and the spin polarizability 7 and 
compare with available information from experiments and from previ- 
ous calculations. Concluding with a critical discussion of our results, 
we point out the necessity of a future 0{e^) calculation. 



* Research supported in part by the National Science and Engineering Re- 
search Council of Canada, by the U.S. National Science Foundation and by 
Schweizerischer Nationalfonds. 



1 Introduction 



Understanding of the implications of QCD within the regime of low energy 
physics has during the past decade become accessible via the technique of 
chiral perturbation theory (ChPT) Initial applications were in the arena 
of Goldstone boson interactions together with a few calculations in the 
baryon sector 0, ^, using relativistic baryon ChPT. In recent years use of 
so-called heavy baryon methods ||^, has generated much interest in calcu- 
lations involving baryons and a great deal of work has been done studying 
strong, weak, and electro-magnetic physics in the near threshold region [0]. 
In this work we will focus on SU(2) heavy baryon chiral perturbation theory 
(HBChPT), which has become the most fully developed sector within baryon 
ChPT 0, I, P]. Thus- far, consistent extension to higher energies in SU(2) 
HBChPT has been limited by treatment of the important A(1232) resonance 
only in terms of its contribution to the various counterterms which arise in 
such calculations. The technique by which to address this deficiency was 
developed some time ago |]lOl by including the A(1232) as an exphcit de- 
gree of freedom in a chiral perturbative scheme^. Recently, a reformulation 
of this formalism which allows for a systematic and explicit calculation of 
higher order terms in an expansion of soft momenta, the pion mass and the 



mass difference niA — rriM has been given, [12, 13| and in this note we apply 



this technique to the problem of forward nucleon Compton scattering and 
the polarizabilities of the nucleon. Nucleon Compton scattering is an area 
of research which has recently received a great deal of attention, both exper- 
imentally and theoretically and in the next section we review the status of 
such work. In section 3 we give a brief introduction to the formalism neces- 
sary to include the A (1232) in chiral calculations, and in section 4 apply this 
to evaluate delta contributions to — 7 scattering, examine its influence on 
the polarizabilities and give a critical discussion of our O(e^) results. Finally, 
in a concluding section 5 we summarize our findings. 



^ Using the formalism of Butler and Savage have given an estimate of the 
contribution of the spin 3/2 resonances to the electric and magnetic polarizabilities of 
the nucleon. However, this calculation was performed in SU(3) and made a number of 
approximations, so that a direct comparison with systematic SU(2) work is not possible. 
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2 Compton Review 

To lowest order the spin-averaged amplitude for Compton scattering on the 
nucleon is given by the Thomson amplitude 

Amp = -^6 -6' (1) 

where Q, M represent the nucleon charge, mass and e, e' and = 
{uj,k),k'fj^ = {uj',k') specify the polarization vectors and four- momenta of 
the initial,final photons respectively. In next order are generated contribu- 
tions arising from electric and magnetic polarizabilities — oe and Pm — which 
measure the response of the nucleon to the application of quasi-static electric 
and magnetic fields 

Amp = e ■ e' (^-j^ + uu' A-na^ + e x k ■ e x k' Att^m + 0{uj^) . (2) 

The associated differential scattering cross section on the proton is given by 

da _ / V/uj'Vll 
dVL ~ \AtiM) \uj ^ 



-(1 + cos^ 



e2 



(3) 



Thus by measurement of the differential Compton scattering cross section 
one can extract the electric and magnetic polarizabilities, provided 

i) the energy is large enough that these terms are significant with respect 
to the Thomson contribution but 

ii) not so large that higher order effects dominate. 

This has been accomplished for the proton in the energy regime 50 MeV< 
uj <100 MeV, yielding 

Q^E = (11.6 ±0.6 ±0.6) X 10" W; Plj = (2.6 T 0.6^0.6) x lO'W (4) 
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(Note: In practice one generally uses the results of unitarity and the validity 
of the forward scattering dispersion relation which yields the Baldin sum rule 

0^1 + 'Pli = TT-J —c^LM = (14.2 ± 0.3) X 10- W (5) 

as a constraint, since the uncertainty associated with the integral over the 
photo-absorption cross section (Ttot('^) is smaller than that associated with 
the polarizability measurements.) 

Since the neutron has no charge, such a Thomson-polarizability inter- 
ference experiment is not possible, so alternative methods must be used. 
The most precise measurement involves a recent n-Pb scattering experiment, 
wherein the dipole moment induced in the moving neutron due to the nu- 
clear charge acts back on the Pb nucleus. This second order process and 
the resulting interaction proportional to the electric polarizability can 
be detected in a transmission experiment via the characteristic linear depen- 
dence on the beam momentum k. This experiment was recently performed 
at ORNL and yielded results flEI] 

al = (12.6 ±1.5 ±2.0) x 10" plj = (3.2^1.5^2.0) x 10"^ (6) 

quite similar to those of the proton. (In this experiment only the electric 
polarizability is measured directly. However, using the unitarity sum rule 
result 

o^I + Bm = ^J — <ot(^) = (15.8 ± 0.5) X 10-W (7) 

/TT^ JO UJ'^ 

the magnetic polarizability can be extracted.) 

On the theoretical side, constituent quark-based approaches which rely 
on the sum rulejlGl 

= TT7 <2^ei[n- RcJj > ±2 2^ (8) 

AM .^^ hn - ho 

are bound to fail since the sum over intermediate states component of Eq. (|^) 
is in almost any reasonable model identical for both neutron or proton, lead- 
ing to a prediction 

«E - «S ^ ^(< rl>-<rl>)= 4.2 x lO'W (9) 
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in conflict with the experimental result that > q;^.[|T8| The solution to 
the problem lies in a proper treatment of the pion degrees of freedom of the 
nucleon, which is suggested by the feature that the leading pion loop contri- 
butions to neutron and proton are identical, in agreement with experiment. 
The problem can best be addressed using the technique of heavy baryon 
chiral perturbation theory 0, within which at 0{p^) one finds a result 

a| = = lOPl, = 10/9]^ = ^/"f = 12.2 x lO^W (10) 

This 0{p^) calculation represents the leading result for ue and Pm in ChPT, 
but gets the qualitative features of the polarizabilities right and even agree- 
ment with experiment ! The results diverge as 1/^7^ in the chiral limit, giving 
support to the idea that at these low energies the photon interacts primarily 
with the long-range pion cloud of the nucleon. In order to understand the 
experimental finding that > the leading ChPT result is not sufficient. 
One must include higher order terms in order to find isospin dependent ef- 
fects and to judge the convergence behaviour of the series^ . A calculation at 
has been performed by Bernard, Kaiser, Schmidt and Meifiner [^]. At 
this order four counterterms are required, which were estimated by BKSM 
by treating higher resonances — including the delta resonance — as very heavy 
with respect to the nucleon. The results of this process are 



a| = (10.5 ± 2.0) X lO^W; (31^ = (3.5 ± 3.6) x lO^W 

= (13.4 ± 1.5) X 10-W; = (7.8 ± 3.6) x lO^W (11) 



where the uncertainty is associated with the counterterm contribution from 
the A and from K, rj loop effects. A very interesting aspect of this C(p^) 
calculation lies in the fact that it identifies a mechanism to counter the large 
positive contribution on Pm due to A(1232) resonance exchange, which is 
a well-known problem in calculations of the magnetic polarizability via ef- 
fective lagrangians BKSM found that at 0{p^) (negative) A^vr-loop 



contributions can essentially balance out the (positive) Delta effects hidden 
in the counterterms ! Nevertheless, we observe that the uncertainties of their 
(9(]9^) calculation are quite dramatic so that real understanding of the ChPT 

^Higher order corrections to the 0{p'^) HBChPT results for the polarizabihties have 
been discussed by L'vov |19||. 
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predictions for the polarizabilities will require more work. In particular, the 
large uncertainty in (3m is mainly related to poorly known couplings involving 
A(1232), which was used to determine some of the counterterms via "reso- 
nance saturation". If one retains A(1232) as an explicit degree of freedom 
in the chiral calculation, as done in our 0{e^) calculation of section 4, one 
can in principle determine all delta couplings of interest from (other) exper- 
iments in a systematic fashion. Furthermore, one does not limit oneself to 
the narrow scope of "resonance saturation" for the spin 3/2 contributions. 
We refer to section 4 for further discussion of this point. 

Finally, with respect to the magnetic polarizability, we note that the 
simple quark model can provide a basic understanding of experiment. The 
prediction! 



J^^. In \ |2 



+ — — - (12) 

involves a substantial diamagnetic recoil contribution 

.dia_/ -10.2 X 10- W p 

~ \ -8.5 X 10~W n ^^"^^ 

which when added to the large paramagnetic contribution due to the 
A(1232)|T1 

.A _/ +12X10-W p 

~ \ +12 X 10-W n ^^^> 

gives results in basic agreement with the experimental findings. Hence, it is 
clear that proper inclusion of the A degrees of freedom is essential. 

When spin-dependence is included the situation becomes somewhat more 
complex. In order to simplify the present exploratory analysis, we restrict our 
attention to forward scattering, in which case the amplitude can be written 
as 

Amp = fi{uj) i ■ e' + ia ■ {i' X i) u f2{uj) (15) 
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where fi{uj),f2{uj) are both even functions under crossing — u — >• —u — and 
are therefore functions only of a;^. In terms of our previous notation we have 



/i(^) = + («E + Pm)u;' + ... (16) 



while in the case of f2{oi!) we can make a similar expansion P2|l 



/2(^) = -|^+7^^ + ... (17) 

where k is the anomalous magnetic moment of the target and the new struc- 
ture 7 is the "spin-polarizability." As in the case of fi{uj) the form of the 
leading term in the expansion is dictated by rigorous low energy theorems, 
while the u"^ correction represents a probe of hadronic structure. One differ- 
ence between the spin-dependent and spin-averaged amplitudes, however, is 
the asymptotic behavior as oo. The better behavior of /2(w) suggested 
by Regge theory allows one to write an unsubtracted dispersion relation in 



this case, leading to the strictures[^ 

vre^K^ duj 



nM2 I L^-(^) 

2M^ Jo uj 



-Jo (18) 



^ An 

where cr±{uj) are the photo-absorption cross sections for parallel and anti- 
parallel alignments of photon and target helicities. The first of these 
relations — the Drell-Hearn-Gerasimov (DHG) sum rule — has received a good 
deal of attention recently. On the experimental side, efforts are being 
mounted to measure the spin-dependent structure function f2{uj) directly, 
thereby confirming the prediction of the low energy theorem. [^] However, 



this has not yet been achieved. On the theoretical end, there have been a 
number of attempts to evaluate the dispersive integral of the DHG using 
what information currently exists for the photo-absorption cross sections. 
The existing data set is incomplete in that hehcity-dependent cross sections 
have not yet been measured. Thus one uses existing multipole decomposi- 
tions from unpolarized experiments in order to perform the analysis. Such 
decompositions are available, however, only in the single pion production 
channel so that above the two pion threshold model-dependent assumptions 
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must be made. With this caveat present results are somewhat higher than 
predicted by the low energy theorem 



°° du. 



(J. 



UJ 



2M2 



multipole 



0.167GeV"2 

0.208GeV-2 

0.210GeV-2 
0.182GeV^2 



Karliner[25] 
Workman/ Arndt 
Burkert /Li 



(19) 



and there have even been arguments made that challenge the assumptions un- 
der which the DHG form was derived. ||2^ However, resolution of these prob- 
lems awaits reliable helicity-dependent cross section measurements, which 
should be available in the near future.^ 

Using the same multipole analysis it is possible to evaluate the dispersion 
integral involving the spin dependent polarizability, yielding! 
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47r2 Jo 



(^)] 



multipole 



-1.3 X lO-^fm-^ 
-0.4 X lO-^fm"^ 



P 

n 



(20) 



For a detailed discussion of the contributions of the various nucleon res- 
onances to the sum rules of Eq.(|18]) and the problems of the constituent 
quark model in describing the DHG sum rule we refer the reader to the 



review article by Drechsel |3C 



These results are in dramatic disagreement with the 0{p^) ChPT 
predictions P] 



7p = 7n 



+ 4.4 X 10-^fm" 



(21) 



Unlike the case of a^; and Pm, at present there exists no corresponding O(p^) 
ChPT calculation that reconciles this discrepancy. Presently the best infor- 
mation about the sub-leading behaviour of 7 comes from a relativistic baryon 
ChPT calculation at the one-loop level f^, yielding 



7;- 



loop 



+2.2 X 10"^ fm^ 



l-loop 

in 



+3.2 X 10"^ fm^ 



(22) 



■^A second area of interest is in the generalization of the DHG form to include the deep 
inelastic region and its connection with the integrated spin dependent structure function 
gi{k^). However, we shall not discuss this issue herein. 
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This calculation does not resolve the discrepancy with the multipole analysis 
results of Eg . (pOD . However, it is known from phenomenological considera- 
tionsQthat A(1232) makes a major contribution of opposite sign. In section 
4.2 we present a systematic chiral calculation of 7 with nucleons, deltas and 
pions as explicit degrees of freedom. 

Having given a brief summary of current research in this area we now 
proceed to outline the formalism which allows inclusion of the A(1232) in a 
consistent chiral power counting framework. 



3 Small Scale Expansion in HBChPT 

The subject of SU(2) Heavy Baryon ChPT of nucleons and pions has been 
well developed in recent years 0- In the conventional form one organizes 
the calculation according to an powercounting, where p denotes ei- 

ther a soft momentum or the pion mass m-n- All nucleon resonances, strange 
particles, vector mesons, etc. are integrated out, i.e. they only appear in 
higher order contact interactions. For near threshold processes this program 
has proved highly successful-for an outstanding review of the field we refer 
the reader to ref. [0. However, phenomenologically it is known that in the 
nucleon sector the first nucleon resonance-A(1232)-lies very close to the nu- 
cleon and can exert its influence even to processes at very low energies. This 
special situation in the baryon sector has early on ||lOl prompted suggestions 
to keep the spin-3 / 2 baryon resonances as explicit degrees of freedom in the 
chiral lagrangian in order to include all relevant physics and to improve the 
convergence of the perturbation series. 

If one retains A(1232) as an explicit degree of freedom in SU(2) HBChPT, 
one is faced with an additional dimensionful parameter A = Ma — Mjv, 
which corresponds to the mass splitting between the nucleon and the delta 
resonance in the chiral limit. Phenomenologically A is a small parameter of 
300 MeV which, unlike remains finite in the chiral limit. As was shown 
in refs. [|rn|, O, O one can nevertheless set up a consistent field-theoretic 



HBChPT formalism provided one organizes the calculation according to an 
0{e^) powercounting, where e now denotes a small scale of either a soft 
momentum, the pion mass or the mass splitting A. Our formalism has been 
constructed in such a way, that an 0{e^) result automatically contains any 



*This has also been pointed out by the authors of 
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result in |^ plus additional terms involving the A(1232) resonance 
which might be present to the order we are calculating. 

In our calculation of Compton scattering below, which is done to 0{e^) 
in the small scale expansion, we shall need only the lowest order vertices, 
the propagator involving the A(1232), as well as possible counterterm con- 
tributions to NNj, NN'j'y, and A^A7. Details of how this can be achieved 



in a general 1/M expansion are given in ref. ||T2|, |T3[. Here we list only the 
minimal results necessary for the present calculation. 

The systematic 1 /M-expansion of the coupled iVA-system starts with the 
most general chiral invariant lagrangian involving relativistic spin 1/2 (iPn) 
and spin 3/2 (ip^) fields ^ 



C = Cn + Ca + (Can + h.c.) (23) 



with 



+y^7.. f'l^ + f (7M^if + u^luH + f 7m ^^^757.] + - 

Cn = i^N [i - Mn + ^ ii^i^N + ■■■ 

^AN = QttNA i'i {9,,u + ^7^*7;^) W^ipN + •■■ (24) 

We have only displayed the leading order terms, the dots denote terms with 
more derivatives or insertions of the light quark mass matrix. It is understood 
that they are included up to the order one is working in the e-expansion. 
Following the conventions of SU(2) HBChPT in the spin 1/2 sector 0, we 
have defined the following structures 



Hii order to take into account the isospin 3/2 property of A(1232) we supply the Rarita- 
Schwinger spinor with an additional isospin index i, subject to the subsidiary condition 
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1 r 



U; 



w, 



T U,. 



—Tr 
2 



U 



d^U - i(v^ + a^)t/ + iU (v^ - 



2 ( ^ ^ ^ 

U = exp —T ■ IT 



(25) 



where v^, denote external vector, axial- vector fields. 

The "light" degrees of freedom of the spin 3/2 field, which are retained 
in the effective low-energy theory, are identified as 



(26) 



where we have introduced a spin 3/2 projection operator for fields with fixed 
velocity 

^^'^ - 9i.u - -lt,lu - - {^It^Vy + ^^m7v ^) ■ (27) 



The remaining components. 



G;(x) = {g,. - P+ Pj/')^J C(a^) exp(zMT; ■ x), 



(28) 



are shown to be "heavy" p2| , |T3| and are integrated out. 

Rewriting the lagrangians of Eq. (|23|) in terms of the spin 3/2 heavy baryon 
components and G^, and the corresponding "hght" and "heavy" spin 1/2 
components A^, h, defined as 



N{x) = P/ ipN exp{iMv ■ x) 
h{x) = P~ ipN exp{iMv ■ x), 



(29) 



the general lagrangians take the form 



Ln 
Lan 
La 



NAnN + (hEnN + /i.e.) - hCNh 

TAanN + GBanN + JiVnaT + JiCnaG + h.c. 

TAaT + (GBaT + h.c) - GCaG. 



(30) 
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Matrices A, B, C admit a low energy scale expansion-^, B start at order e^, 
whereas C has a leading term of order 1. This allows to perform a systematic 
1/M-expansion, following an approach developed by Mannel et al. in the field 
of heavy quark physics , which was later applied to spin 1/2 HBChPT 
by Bernard et al. ^j. The result of this procedure is the effective action for 



the coupled A^A-system [0, 



S^s = J d^x {tAaT + NAnN + [tAanN + h.c] } 



(31) 



with 



An 
Aan 



and 



^A + Io'DnaIoCn^'^na + 7o'Ba7oCa^'Ba 
An + 'JoB^nIo^n^^n + IoB^anIqCa^Ban 
Aan + 'Jo'DnaIoCn^Bn + IoB^^'JoC^^Ban 

Cn = Cn — Cn aC a^IqCn aIo 
Bn = Bn + CnaCaBan 
Vna = T^NA + CnaCa^Ba- 



(32) 



(33) 



The vertices relevant for our calculation can be read off directly from 
Eg ■ (|3ll) . Further analysis below will be divided into two classes of contri- 
butions, one- loop graphs and Born graphs. As one can see from Eg. (|3T|) , 
any calculation at 0{e^) which has a nucleon in the initial (and final) state, 
completely contains a 0{p^) calculation in the formalism of ||^. Furthermore, 
it shows unambiguously, which additional diagrams involving A(1232) must 
be appended. 



3.1 Loop Graphs 



To order e^, only one- loop graphs with vertices of order e must be considered. 
Thus, the pieces in Eg.(^l]) we need are An\ A^^ and Aa^n- More explicitly, 
these are found as 



•^AN 

A^^ 



iv ■ D + qa S ■ u 

g^NAW^ 

-'iv D'^ - A 6'^ +giS- u'^ 



(34) 
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where S*^ denotes the Pauh-Lubanski spin vector [^. From matrices W'^^ 
and A^l^ we determine the SU(2) HBChPT propagators in momentum space 
with soft momentum^ r^ = p^ — Mv^ : 

S^'\vr) = 

V ■ r + IT] 

«^("-^' - ..;!^./ t3/.. (35) 

with P^l^ denoting the spin 3/2 HBChPT projector in d-dimensions ||13| 

Pll^ = 9,u - V,, + ^ S,S, (36) 

and 

el3/2 = - \ r'r^ (37) 

being the corresponding isospin 3/2 projector. From Eq.(^) one can see 
that the deha propagator counts as a quantity of order in our expansion 
scheme. 

At present, we have no systematic determination of the couphng constant 
Q-kNIi within the small scale expansion. For the time being we rely on the phe- 
nomenological analysis of Davidson, Mukhopadhyay and Wittman (DMW) 
m, yielding 

^.Ta^ = - ^?.Ta^ - 1.5 ±0.2. (38) 

Finally, choosing the velocity vector = (1, 0, 0, 0) and working in the 
Coulomb gauge v ■ e = v ■ e' = 0, we conclude that we have to calculate 18 
loop diagrams, displayed in Fig. 3 and Fig. 4. Details of the calculation are 
given in Appendices B and C, the results will be discussed in section 4. 



^Working to 0{e^) and in the absence of any 0(e) Born diagrams due to our choice 
of gauge, we always can identify the large mass M with the physical mass of the nucleon 
Mat. Mass renormalization only becomes important in an ©(e*) HBChPT calculation of 
Compton scattering. 
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vertex 



lagrangian 



N 



0(6) 

0(e) 



-fNN 
-fNA 



0(e2 
0(e2 



-fNN 
-fAN 
e" ) 77A^A^ 



A 



in Coulomb gauge 



and 7oi3iJ^ So4°^ 



7oi3_ 



(1) /-(O) 



-1^(1) 



^ »(2) ^(0) 



N 



h.c. 



Table 1: Vertices for Born graphs in Forward Compton Scattering 



3.2 Born graphs 

The Born graphs contributing at 0{e^) Compton scattering are shown in 
Figs. 1 and 2, with the pertinent vertices given in Table 1. The structures 
not involving the delta resonance can be taken from ref. 0^. We note that 
our use of the Coulomb gauge dramatically reduces the number of diagrams. 
Also, the possible diagram of Fig. 2. Id with an anomalous vr" 77 vertex 
does not contribute in forward direction. Due to the fact that the photo- 
excitations of A(1232) begin with the Ml transition, there exists no 7A^A 
vertex at 0{e). Consequently, there is no 1/M-corrected interaction of this 
type at O(e^) in Table 1. However, the relativistic counterterm lagrangian 

mm 

ibi ^„ , .1 



C 



NA 
c.t. 



(39) 



provides the Ml 7A^A transition strength and leads to an O(e^) structure: 



AN 



with the chiral field tensor being defined as 



(40) 



± 



U 



± u u 



^In general this is not the case. As will be shown in (l3| the "heavy components" of 
the relativistic spin 3/2 field modify the counterterms of the NN lagrangian starting at 
0{€^). In our specific case the NN vertices are unchanged. 
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pL,R pL,R 



F« = v^ + a^, F^ = v^-a^. (41) 



We note that the off-shell parameter y in Eq. ( |5UD does not contribute at order 
when going to the effective heavy baryon lagrangian-it will only enter at 
0{e^). However, since the Delta propagator Eq.(p5D counts as order e~^, 
and we effectively have no 0{e) vertices around, we do not have to consider 
g3 ^]Sf]Sf or 7A^A vertices to the order we are working. Furthermore, the 
absence of an 0{e) '-/NA vertex is also responsible for the fact that the 0{e^) 
two-photon seagull term does not get renormalized by Delta interactions. 
The only Born diagrams involving A(1232) are therefore s- and u-channel 

(2) 

resonance exchange with vertices from A)^^ (diagrams 2.2a,b in Figure 2). At 
present, the magnitude of the finite 0{e^) counterterm hi is not known very 
accurately. Until one has identified a suitable hi dependent observable which 
has both been calculated in the small scale expansion and been reasonably 
well measured!, we will employ a phenomenological relation found in ref. 



321, which in our convention reads 



TT) 

2^ ^ - ( 2-5 ± 0.35 ). (42) 

For the numerical estimate we have used g^^A Eq. (p8|) . This completes 
the background necessary for the present application. 



4 Forward Compton Scattering and Nucleon 
Polarizabilities 

Having outlined our formalism, in this section we investigate the influence of 
the A(1232) resonance on the nucleon polarizabilities. We restrict ourselves 
to the case of forward scattering, which provides information on the electric 
polarizability cue, the magnetic polarizability (3m and the spin polarizability 
7. As mentioned above, for the spin-averaged quantities aE and Pm the 
results involving only nucleon and pion degrees of freedom are already known 
to 0(p^) ||20| in the chiral expansion. Here we present a systematic analysis of 

®Work along these lines is under way; V. Bernard, T.R. Hemmert, J. Kambor and 
U.-G. Meifiner, in preparation. 
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the A(1232) contributions to aE, Pm and 7 at 0{e^). A complete calculation 
of all polarizabilities at O^e"^) including pion, nucleon and Delta degrees of 
freedom will be the subject of future work. 



4.1 Spin-averaged Forward Compton Scattering 

Working in the gauge — e ■ v = — the spin-averaged Compton tensor in for- 
ward direction Q^^, can be written as 

,1. 



[ e' • e U{uj) + e' ■ k e ■ k V{uj)] , (43) 



where k is the four-momentum of a photon with energy u = k-v. e (e') refer to 
the polarization vector of the incoming (outgoing) forward scattering photon 
and T^u{v, k) is the Fourier-transformed matrix element of two time-ordered 
electromagnetic currents 



T^.iv,k) = j d'x e*^(iV(t;)|r[ J7(x)jr(^) ] \N{v)) (44) 

In the spin-averaged case all the information about the low-energy structure 
of the photon is contained in just two function^ Uiy^), Viyj)- However, there 
exists a structure- independent constraint with respect to U{ui), stating 
that in the limit of zero photon energy one has to obtain the Thomson result 
of Eq.(g) 

[/(O) = , (45) 

where Z refers to the charge number of the Compton target and M is its 
mass. Furthermore, U{uj) has to be even under crossing symmetry^, i.e. 
U{uj) = U{-uj). 

Keeping these two non-trivial constraints in mind we split up the calcu- 
lation of U{uj), V{uj) into four separate components: 

U{uj) = f/^(cj)^°™ + f/jv(cu)'°"P^ + [/A(cj)^°'^" + [/A(cu)'°°P^ 

V{uj) = V]v(cu)^°'-" + V]v(u;)'°°^' + Va(cu)^°™ + Va(^)'°°^' (46) 



^The auxiliary function V{uj) can be eliminated for real photons. Nevertheless one can 
obtain information about the magnetic polarizability /3 from it. 

-"[/(w) is related to the function of Eq.(p^) via fi{w) = U{lu). 



lOj 
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We start with the calculation of the nucleoli Born contributions to 
U{uj), V{u!). In the Coulomb gauge e ■ v = there exist non-zero contri- 
butions at O(e^) and 0{e^), as shown by the diagrams in Figures 1 and 2. 
To 0(e3) we find 

with Un{^^) solely stemming from the 0{e^) seagull diagram of Figure 1 
and Vn{uj) arising from diagrams 2.1a and 2.1b of Figure 2. We also note 
that ?7iv(0)^°™ satisfies the Thomson limit Eq.(^5|) for proton {Z = 1) and 
neutron {Z = 0) targets, as expected. We therefore conclude that 

UNiOY""^' = t/A(0)^°"" = f/A(0)'°°P" = , (48) 

and this will serve as a powerful constraint and check on our calculation. 

The C(e^) nucleon loop contributions to the spin- averaged functions can 
be obtained from the nine diagrams shown in Fig. 3 yielding]^, as detailed 
in Appendix B 

Note that the loop effects are isospin-independent {i.e. identical for neutron 
and proton) to this order. We have also checked that U]\f{uj = 0)'°"^'* = 
(Eg. (^8])), by carefully analyzing the dimensionality dependence of the nine 
loop amplitudes as given in Appendix B. 

Next, we evaluate A(1232) Born contributions to U{uj), V{uj). At 0{e^) 
we find two contributing diagrams, as shown in Fig. 2. They yield 

TT /, ABorn 8^1^^ ^ 



9M^ A2 - 



y-i-r'" = - JftiA-. (50) 



11 



Eq.(M9) agrees with the 0{p^) resuh of 
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where bi is the Ml AA^7 couphng of Eq.(^). Again, the Thomson constraint 
[/a(0)^°™ = holds, as required by Eq.®. 

Finally, the corresponding delta loop contributions can be found from the 
nine diagrams shown in Figure 4. Analyzing the invariant amplitudes given 
in Appendix C, we find 



I loop 



547r2F2 



9A 



9m^ 



A2 



(A2 - m2)3/2 



logi? 



A2 



logi? 



mt 



.loop 



UttNA 



547r2F2 



logR+ 0{uj^) 



mi 



with R defined as 



R 



+ 



\ 



A2 



mt 



1. 



(51) 



(52) 



Again we note the validity of the Thomson stricture U^i^^ = 0)'°°^^ = 0. 

In order to extract the nucleon electric and magnetic polarizabilities, we 
now define the nucleon Born term subtracted quantities U{uj), V^u) 



IJ{UJ) 

and make the connections 



U{uj) - f/^(a;)^°™ 
Viuj) - Viv(^)^°™ 



^^f/(-)Uo 



(53) 
(54) 



(55) 
(56) 
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Adding up all three contribution^ one finds aE and I3m to 0{t 



CtE 



+ — 



59l 



1 



4:71 967rF2 m, 




+ 



6 9nNA ^ 



An 547rF2 



TT 



9A 



A2 



+ 



A2 - lOm^ 
(A2 - m2)3/2 



logi? 



(57) 



[12.2 (N-loop) + (delta-pole) + 8.6 (delta-loop)] x 10"* W 



+ 



+ 



bl 1 



47r 9 A 

<2 yTrAfA 



47r 547rF2 /^^ 



2 vr 



logi? 



(58) 



= [1.2 (N-loop) + 12 (delta-pole) + 1.5 (delta-loop)] x 10"* fm^ 

In assessing these results, we observe that from the delta pole terms there 
exists a significant contribution to the magnetic polarizability but none to the 
corresponding electric polarizability at this order ! The strong effect on is 
not surprising, as the large Ml nucleon-delta coupling is known to contribute 
substantially to the magnetic polarizability [p!7| , pl[] . Also, we note that at 
0{e^) the Avr-loop contributions to (3m are of the same size as the A^vr-loop 
effects, which is unexpected. Furthermore, and perhaps most surprisingly, 
we find a large Avr-loop component in a^; ! Even though the numerical values 
of the delta contributions in Eqs.( |^ , |55D are understood to have sizable error 
bars due to quadratic dependence on the presently poorly known couplings 
Q-kNA and bl, it is clear that one cannot expect to achieve agreement with 
experiment at this order of the calculation-the effects of A (1232) are sizable 
and strongly renormalize the A/^vr-loop results. In the 0(p*) calculation of 
[ pO[[ it has been shown that the large delta pole-contribution in Pm cancels 



12r 



2 The A^TT-loop parts of Eqs.(|^,|8|) agree with_the 0{p^) calculation of [|. In SU(3) 
an estimate of spin 3/2 contributions to a^; and Pm was given in [ pi] ]. aE and (3m have 
also been calculated in SU(2) relativistic baryon ChPT [p3[. 
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to a large extent against A^vr-loop effects at that order. In addition to the 
simple delta poles the counterterms of BKSM have been determined via 
"resonance saturation". The large Avr-loop effects that we find at 0{e^) are 
therefore not accounted for[^ in the counterterms of ! Unless one finds 
a cancelation through Avr-loops in a future C(e'^) calculation, one would 
have to conclude that "resonance saturation" via simple pole-graphs is highly 
suspect for counterterms in the baryon sector]^. Future work will address 
this important question. Finally, we note that the results for ctE and Pm are 
isospin independent at C(e^). The surprising experimental results > a^, 
(3^ > /^E cannot be addressed at this order and also warrant an investigation 
at sub-leading order. 

In order to obtain an estimate of the convergence of the perturbation 
series and to check the field-theoretic consistency of our calculation, it is 
useful to perform a chiral expansion of our results for the polarizabilities, 
yielding 



(Xe 



2 _ r 1 /2A 

1 + - log — 

9 Knij, 



9wNA "^TT 



M 




(59) 



97r \mn/ 



(60) 



We note that the long-range pion cloud, which scales as l/m^,, provides 
the dominant singularity in the chiral limit, as expected. Also, one can 
see that decoupling of the delta resonance in the chiral limit is manifest. 
Finally, we note that the leading order p^-terms in the chiral expansion of 
the Avr-loop contributions is already a good approximation to the full order 
expressions. The contributions to a^; and Pm are thereby changed by 15 
and 7 % respectively, indicating that the bulk of the effects due to A(1232) 

^■^In particular, the counterterm contribution to has been estimated to be Sa w 
2.0 X lO-'^ W [|0) ! 

^"^Of course these counterterms can still be determined from other experiments. A 
breakdown of the "resonance saturation" hypothesis only means that the numerical value 
of the counterterms cannot be understood within a simple resonance model. 
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are obtained at order p*^ in the chiral expansion ! Of course we emphasize 
that these considerations are only an indication of what might happen at the 
next order, future work will address these issues. We now move on to discuss 
the case of spin-dependent quantities in Compton scattering. 



4.2 Spin-dependent Forward Compton Scattering 

In the presence of spin-dependence the Compton tensor for forward scattering 
of real photons Eq. (|43|) from a spin 1/2 particle has to be expanded by extra 
(spin-dependent) structures. Choosing = (1, 0, 0, 0) and again working in 
Coulomb gauge e ■ f = e'- = we can write 



- i' -e U{uj) + e -ke-k V{u) 
+i 00 W'^^\uj) a-{e' xe)+iuj W'^^\uj) e-ka-(e' xk) 

W^'^^{uj)(^ -ko ■(exk)\ , (61) 

where uj = v ■ k denotes the energy of the forward scattering photon with 
four- momentum /c^. In the following we enforce the transversality condition 
e ■ k = e' ■ k = and therefore only work with the auxiliary function]^ 
W^^\u), which contains the information about the spin 1/2 structure of the 
target nucleon. Having factored out an extra u, we note that W^^\uj) is even 
under crossing (i.e. receives contributions from 

four different sources, analogously to Eq. (|46|) : 

W^'\uj) = W4'^(cu)^°™ + W^^\ujy°°P' + Vri'V)^"™ + W^^\ujy""P' (62) 

First we calculate the nucleon Born contributions to W^^\uj). They arise 
at at 0{e^) from the four diagrams 2.1a-d in Figure 2. One finds 

W^^\ur- = {4 (1 + r,) + ^1(1- n)} , (63) 

where („) corresponds to the anomalous magnetic moment for a proton 
(neutron) target. As noted before, the contribution from the anomalous 

i5pj.Q]-^ pi/(2)(^) and W^^^Lo) one can obtain the third-order spin polarizabilities 73, 74 
of ref. p^ . However, we relegate this analysis to future work where we will study the 
complete set of third-order spin polarizabilities of the nucleon. 
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process 7r° — > 77 (diagram 2. Id in Figure 2) vanishes in the forward direc- 
tion. However, the sum of the other three diagrams satisfies the famous LET 
due to Gell-Mann, Goldberger and Low ||22|, with the two-photon seagull 
diagram 2.1c arising from 1/M^ corrections (see Table 1) playing a pivotal 



role. Relating 



(1), 

N 



[Born 



of Eq.(|63|) to the function /2(u;) of Eg. ([T5|) via 



/2(ct;) = e^/47r W^^\lj), one finds in the limit of zero photon energy 



/2(0) 



8iT Ml 



(64) 



Having satisfied the low energy constraint we conclude that 



\Pole 



1^«(0)'°°P = 



(65) 



which constitutes a non-trivial check on our calculations. 

The 0{e^) A^vr-loop contribution to W'^^\uj) can be obtained from just 
six diagrams (3a-3f) of Figure 3. With details of the calculation given in 
Appendix B we findQ 



(1), 



N 



UJ 



\ loops 



9\ 



2FI Vl-n'^ml 



+ 0{uj' 



(66) 



so that the LET constraint W^^\^)^°°'p = is satisfied. 

The Oie"^ delta Born contributions to the spin-flip function ^^'^^•'(a;) are 
given by diagrams 2.2a and 2.2b of Figure 2, yielding 



W 



(1), 



UJ 



UJ 



.Born _ _ ^"1 

A2 - u;2 



(67) 



Again, the LET stricture Eg . (|65|) holds explicitly. 

Finally, the O(e^) Avr-loop contributions to W'^^^iyj) arise from diagrams 
4a-f in Figure 4. With details of the calculation given in Appendix C, we 
find 



W^^[uj)^°°^ 



91na 
547r2 F2 



3mlA 



(A^ 



m 



2)5/2 



\ogR 



A^ + 2ml 
(A2 - ml)' 
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Eq.(pq) agrees with the 0{p^) result of | 



(68) 
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where R has been defined in Eg. ([5^) . The d-dependence of the pertinent 
amphtudes in Appendix C-induced by the spin 3/2 propagator of Eq.(^)- 
leads to a highly complex cancelation pattern among the diagrams, ultimately 
yielding wjl\oy°°P = as required by Eq.(|^). 

In analogy to the spin- independent discussion in section 4.1 we introduce 
the nucleon Born-term subtracted function 

(1) (uj) = W^^^ (uj) - Wi^^ (69) 

and obtain the spin polarizability 7 of the nucleon 



7 



Stt dtu 



W^'\co)l=o, (70) 



which was defined via Eq.(^). 

Summing the three contribution^ contained in W^-^^lu) one finds 7 to 



o2 „2 



Att 247r2F2 ml 

ynNA 



+ 2ml A 



(A2 _ miy (A2 - m2)5/2 



logi? 



(71) 



= (4.5 (N-loop) - 4.0 (delta-pole) - 0.4(deha-loop)) x 10~^ fm^ . 

We note that as in the case of a^; and Pm there is no isospin dependence at 
0{e^). Also, Ayr-loops are only playing a minor role in the spin polarizability 
at this order. Finally, the large positive contribution of the iVyr-loops is 
nearly completely canceled by the delta Born graphs. In the case of the 
spin polarizability we cannot yet make a direct experimental comparison, as 
explained in section 2. However, comparing with the sum rule value given 
in Eq.(|T8|) we see that the A(1232) contribution goes in the right direction 
but is not large enough in (negative) magnitude in order to bring about 
experimental agreement. We also note that at present we do not have a 
(9(p^) calculation in HBChPT with which to compare. However, the one-loop 



17 



Our findings agree with the Nn- and Att-Ioop results given in [| 
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relativistic baryon results of Eq.(p2D already give an indication that even if 
one extends the 0{p^) HBChPT result of Eq.(^T]) to the next order, one will 
not be able to describe 7 without keeping A(1232) as an explicit degree of 
freedom in the theory. Studying the chiral limit of the spin polarizability 
supports this viewpoint. We find 



Att 2167r2 F3 ml 



9gl - 6bl 



(47rF.) 



A2 



(72) 



One observes that the 1 / singularity due to the pion cloud of the nucleon 
remains leading in the chiral limit and that delta decoupling holds. We also 
note that the Avr-loop contributions increase in magnitude to —0.6 10~^ fm^ 
but remain small. Most importantly, Eq.(|72D shows us that all delta effects 
only start contributing at 0{p^) (!) in standard SU(2) HBChPT. The spin 
dependent polarizability of the nucleon is therefore a prime example of how 
a theory with explicit delta degrees of freedom can dramatically improve 
the convergence of the perturbation series. In light of the cancelations in 
Eq. ([7T|) we are encouraged that at C(e'') one may be able to achieve a good 
understanding of the spin polarizability with the help of the small scale ex- 
pansionQ. 



4.3 Critical Discussion of 0{e^) Results 

Before concluding it is useful to discuss some of the arguments of which we are 
aware why an 0{e^) calculation can perform poorly on the phenomenological 
side and on the other hand show all the correct field-theoretic constraints of 
LETs, heavy mass decoupling, correct chiral limit etc. 

i) The heavy baryon propagator for nucleons of Eq.(|35D receives "propaga- 
tor corrections" via 1/M vertices in matrix A^^^ of Eq. (|32D , which in one- loop 
diagrams only begin to show up at 0{e^). From a relativistic viewpoint these 
1/M corrections correspond to effects pertaining to the lower components of 
a Dirac spinor. 

There certainly have been many attempts in the hterature to cah the delta to the 
rescue in order to get reasonable numbers for the sum rules of Eq.(]l8|). We emphasize again 
that our results in Eq.([7l|) follow from the systematic small scale expansion in HBChPT, 
as laid out in section 3. 
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ii) For the heavy baryon propagator of a spin 3/2 particle there are anal- 
ogous 1/M induced propagator corrections in matrix A^^^ of Eg. (|32D . In this 
case they have the additional effect of correcting the pure spin 3/2 projec- 
tor of Eq.(|53|) by bringing in some information about the off-shell spin 1/2 
components of a relativistic spin 3/2 particle. Phenomenologically it is well 
known that the spin 1/2 components of the relativistic spin 3/2 propagator 
can play an important effect in some observables |3^. Again, at the one- loop 
level these effects only start showing up at 0{e^). 

iii) In 0{e^) calculations the incoming (and outgoing) soft momentum of 
the baryon is usually a higher order effect. Only at the next order does this 
momentum fully contribute and the size of these "recoil corrections" can be 
quite large for some observables. 

iv) Isospin is usually not broken in calculations at the 0{e^) level which 
can be a problem in some cases. 

v) Often one is calculating observables which are dominated by loop ef- 
fects. Therefore one is only calculating the leading order of these observ- 
ables, although technically one might have to work to third or fourth order 
in HBChPT. As experience from many calculations in the meson sector of 
ChPT teaches us, one should always consider the leading and the first (at 
least) non- leading order in an observable before one can try to judge the 
quality of a ChPT result. 

This list is certainly not complete. The concerns raised are valid be- 
yond our particular case of the 0{e^) calculation of the polarizabilities of 
the nucleon. In fact, most of these points also apply to SU(2) HBChPT 
without explicit delta degrees of freedom and "leading log" calculations in 
SU(3) HBChPT. For our particular case of interest — the polarizabilities of 
the nucleon — all of these points strongly suggest moving onto the 0{e^) cal- 
culation. 



5 Conclusions 

The Compton scattering process offers the opportunity to probe nucleon 
structure in a relatively clean fashion, via measurement of various "polar- 
izabilities" which probe the nucleon response to quasi-static excitations. In 
particular recent years have seen the measurement of electric and magnetic 
polarizabilities for both neutron and proton and spin polarizability measure- 
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merits should be available by the millennium. On the theoretical side, 0{p^) 
predictions for cxe and Pm within heavy baryon chiral perturbation theory 
give a surprisingly good picture of the experimental situation. However, cor- 
responding predictions for the spin polarizability are not in good agreement 
with the values given from DHG sum rule arguments. When extended to 
0{p^) reasonable agreement is obtained for dE and f^M, but at the cost of 
considerable uncertainty associated with effects such as the A(1232), which is 
included only as a very heavy particle contributing to various counterterms. 
In this paper we have removed this obstacle by treating the A(1232) as a 
specific degree of freedom within the heavy baryon method. Our calculation 
was performed to 0{e^) where e is taken as a soft momentum, as m^r or as 
the nucleon-delta mass difference. We find, perhaps not surprisingly, that 
inclusion of delta effects makes large changes to all three polarizabihties. 

The good agreement with experiment at 0{p^) for a^; and Pm is destroyed 
at 0{e^), whereas the spin polarizability 7 improves dramatically when com- 
pared with currently available sum rule information. Furthermore, we have 
discussed consequences of our HBChPT results in the light of existing 0{p'^) 
and relativistic one-loop calculations. Indeed, we regard our calculation as 
merely preliminary and look forward to extending this work to higher order. 
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A Loop Functions 

We express the invariant amplitudes of Feynman diagrams containing pion- 
loops around a nucleon in terms of J-functions, defined via 



{1 



5 ^IJ,^ui ^fi^u^a^/B 



} 



i J {27rY {v-e-uj- ie){ml -P- ie) 



Employing the following identities 

1 



{Joiuj,ml), 



(73) 



+ perm.)J6(^^,"^^) + •••} 



J2 (a;,my 



d- 1 
1 

d + 1 



(m^ - u;^) Jo rnl^ - to 
(ml - u^) J2 {uj, ml) 



rniuj 



— A 

d ' 



(74) 
(75) 



one concludes that all loop-integrals can be expressed via the basis-function 
Jq. For the A^vr loop-integrals we use |^ 



-UJ 



V - ^'arccos— , (76) 

whereas for the Ayr loop-integrals we employ the analytically continued func- 
tion 



Jo (uj,ml) 



UJ 



m^ 



ALoj + — 1 - 21n— 
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87r2 



/i / 47r2 



^2 _ ^2 



X log 
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In Eqs. ([7^), ([75|) , (|76| ) and ( ffTl) we have used the same conventions as 



(77) 



A, 



2ml \L + 



1 m^ 
-m- 



+ O (rf - 4) 



(78) 
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where we introduced the Euler-Mascharoni constant, 7^; = 0.557215, and the 
scale n in the dimensional regularization scheme we use for the evaluation of 
the integrals. 

Finally, with and J" we define the first and second partial derivative 
with respect to m^, 




(80) 



(79) 
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B Nil Loops in Forward Compton Scattering 



Using the J-function formalism defined in Appendix A and the lagrangians 
of Eg . (13^ one can get exact solutions for the nine iVvr-loop diagrams of 
Figure 3. With (e^) we denote the polarization four-vector of the incoming 
(outgoing) photon with constant four-momentum and energy uo. We find 



Amp^^^ = Cu{r){ 



Ntt 



e ■ e 



Jo{uj,ml) + Jo(-u;,m^)] 



(81) 



+2 [S^,S,] e'^^e^ [Muj,ml) - Jo(-^,<)] } ^(r) 



Amp^l = CM(r){ +4:e-e' dx[j2{ujx,ml) + J2{-ujx,ml)] (82) 
—8 [S^,S,y] t ^t'^ I dx J2{ujx,ml.) — J2{—ujx,m1) 



—2e-ke'-k dx x{l — 2x) Jq{ujx, m^) + Jq(— tux, m^) 

JO L 

+2 [S^, (e' ■ k kf'e'' + e ■ k e' "k") f dx x{l - 2x) 

Jo 

X [jo(u;x,m^) - Jo(-u;a;, m^)] } u{r) 
C u{r)-[ -Ae-e' (d+l) £ dx (l-x) 



J'(!{ujx,ml) + JQ{-Ujx,ml) 



(83) 



+4 e ■ e oj i dx {1 — x)x J^^iux, m^) + (— cux, m^) 

Jo '- 

- e-ke' -k dx {I - x) \%x{2x - 1) + {2x - lf{d - 1) 



X \ j'^{ujx,ml) + j^^-ujx.ml] 

+ e-ke' -kuj^ dx (1 - x) x^ (2x - 1)^ 
Jo 

X [jo{^x,ml) + JQ{-Ujx,ml)'^ } u{r) 
Cu{r){ +2e-e' (d-l) 4{0,ml)} u{r) , (84) 
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with the common factor 



C = i 



2F2 



In order to evaluate the polarizabihties, we expand the nine amphtudes into 
a power series in u and only keep the terms of interest: 
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(88) 



In section 4 we construct the auxiliary functions U{uj), V{ijj) and Ty(u;) from 
Eqs.(|5l]8|). 
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C Att Loops in Forward Compton Scattering 

Using the same conventions as in Appendix B with a new overall factor 



3F2 



we find the invariant amplitudes for the nine diagrams of Figure 4: 
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In analogy to Appendix B we expand the amplitudes into a power-series 
in the photon energy uj. With the definition 
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we find: 
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FIGURE CAPTIONS: 

Wiggly lines to the right (left) denote incoming (outgoing) photons, dotted 
lines denote pions and solid lines represent nucleons. 

Figure 1: O(e^) Born graph in forward Compton scattering. 

Figure 2: O(e^) Born contributions in forward Compton scattering. 

Figure 3: O(e^) A^tt Loop diagrams in forward Compton scattering. 

Figure 4: O(e^) Att Loop diagrams in forward Compton scattering. 
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gure 1: O(e^) Born graph in forward Compton Scattering. 
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Figure 2: 0(e"^) Born graphs in forward Compton Scattering. 
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gure 3: 0(e"^) Ntt Loop Diagrams in Forward Compton Scattering. 
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Figure 4: 0(e"^) Att Loop Diagrams in Forward Compton Scattering. 



